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We study QCD with 2 and 3 flavours of colour-sextet quarks. The 2-flavour theory is a candi- 
date Walking Technicolor theory. Since we are attempting to distinguish whether this theory is 
walking or conformal, we also study the 3-flavour theory, which is believed to be conformal, for 
comparison. We simulate lattice QCD with 2 and 3 flavours of colour-sextet staggered quarks 
at finite temperatures to determine the scales of confinement and chiral-symmetry breaking from 
the positions of the deconfinement and chiral-symmetry restoration transitions. Unlike the case 
with fundamental quarks, these transitions are far apart. For 2 flavours the values of /3 = 6 /g 2 for 



both transitions increase as Ta is decreased from j to g to g, as expected for a theory whose cou- 
pling runs to smaller values as the lattice spacing is decreased. However, for the chiral transition, 
the increase in j3 between Ta = | and Ta = -L is much larger than the increase between Ta = J 

4 we are at strong coupling where 



and Ta = i. This suggests that between Ta = | and Ta 
the theory is effectively quenched, while between Ta = g and Ta = | we are emerging into the 
weak coupling regime. It will require even smaller Ta values to determine whether the running 
of the chiral-transition coupling is controlled by asymptotic freedom and the theory walks, or if it 
reaches a non-zero limit when the transition becomes a bulk transition and the theory is conformal. 
The 3 flavour case at Ta = \ and Ta — I behaves similarly to the 2 flavour case. Since this the- 
ory is expected to be conformal, the interpretation that we are seeing strong-coupling behaviour, 
inaccessible from the weak-coupling limit (continuum) is the most likely interpretation. 
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1. Introduction 

Technicolor theories are QCD-like theories where the techni-pions play the role of the Higgs 
field, giving masses to the W and Z [||, §]. If the choice of gauge fields and fermions is such that 
the coupling constant evolves very slowly - walks - phenomenological problems can be avoided 
when this Technicolor theory is extended to also give masses to the quarks and leptons []3|, f|, ||, gj. 

QCD with < Nf < 3jq flavours of massless colour-sextet quarks is expected to be either a 
walking or conformal field theory. (For a good summary of what is known from perturbation theory 
about SU(N) gauge theories with Nf fermions in various representations, see [j7j|.) For Nf = 3, 
conformal behaviour is expected. The Nf = 2 theory could, a priori, exhibit either behaviour. We 
simulate lattice QCD with Nf = 2 colour-sextet quarks to determine which of these two options it 
chooses. We also simulate the Nf = 3 theory for comparison. 

Our earlier studies of the Nf = 2 theory on lattices with N t =4,6 indicated that the well- 
separated deconfinement and chiral-symmetry restoration transitions move to weaker couplings as 
N t increases from 4 to 6 as expected of the finite temperature transitions of an asymptotically free 
theory [H]. This seemed to favour the 'walking' scenario. We are now extending these studies to 
./V, = 8. Again both transitions move to weaker couplings than those for N t = 6. However, for the 
chiral transition, the decrease from N t = 6 to N t = 8 is far less than that from N t = 4 to N t = 6. 
This raises the possibility that the coupling will tend to a finite limit as N t — > °°. If so, the transition 
is a bulk transition, and the theory is conformal. Otherwise, if this critical coupling vanishes in 
this limit, the transition remains a finite temperature transition and the theory walks. It will be 
necessary to simulate even larger N t to distinguish between these two possibilities. 

We are also performing finite temperature simulations of the N/ = 3 theory on N t = 4, 6 lattices. 
The behaviour of this theory appears very similar to that of the Nf = 2 theory; the deconfinement 
and chiral-symmetry restoration transitions are widely separated, and move to significantly weaker 
couplings as A^ increases from 4 to 6. The main difference is that all transitions occur at stronger 
couplings than for Nf = 2, as expected. Since this theory is expected to be conformal, this suggests 
that the N t = 4,6 chiral transitions are in the strong-coupled domain, which is separated from the 
weakly-coupled conformal regime by a bulk transition, and is hence inaccessible. 

2. Nf = 2 simulations and results 

We are simulating lattice QCD with the Wilson plaquette action and 2 flavours of (unimproved) 
staggered colour-sextet quarks using the RHMC algorithm. We study thermodynamics by simulat- 
ing on a lattice with temporal extent N t = 1/ (7a) where T is the temperature and a is the lattice 
spacing. The spatial extent of the lattice A^ S> A; (We are interested in the limit N s — > °°). 

Our earlier work used lattices with A 7 , = 4 and 6. We are now simulating on lattices with A^ = 8. 
The results we present here are on 16 3 x 8 lattices, with quark masses m q = 0.005, m q = 0.01 and 
m q = 0.02. We perform simulations with /3 = 6/g 2 in the range 5.5 < j8 < 7.4. Away from the 
deconfinement and negative-to-complex Wilson-line transitions, we perform runs of 10,000 length- 
1 trajectories for each (j3,m). Close to these transitions we increase this to 50,000 trajectories. To 
check for finite size effects we are also performing simulations with m q = 0.005 at large /3s on 
24 3 x 8 lattices. 
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Just above the deconfinement transition, the Wilson line shows a strong 3-state signal, where it 
preferentially orients itself in a direction close to that of one of the cube roots of unity. We therefore 
bin our 'data' into 3 bins according to whether the argument of the Wilson line is closest to 0, 27i/3 
or — 2n/3, up until the complex Wilson-line states disorder to a negative Wilson-line state. 

Figure [l] shows the Wilson Lines and chiral condensates as functions of j3 for each of 

the 3 quark masses, for our Nf = 2 runs on a 16 3 x 8 lattice. This 'data' is for the real positive 
Wilson-Line bin only. The deconfinement transition is identified as the point at which the mag- 
nitude of the Wilson Line jumps from near zero to an appreciably larger value. From histograms 
of the magnitudes of the Wilson lines, we estimate that the deconfinement j8, fa = 5.665(10) for 
m = 0.02 and fa = 5.660(10) for m = 0.01. 
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Figure 1: Wilson Lines and chiral condensates for 
the positive Wilson Line states, on a 16 3 x 8 lattice. Figure 2: Chiral susceptibilities on a 16 3 x 8 lattice. 



In figure ||] we plot the disconnected chiral susceptibilities for each of the 3 quark masses. 
For both m = 0.01 and m = 0.005, the peaks are at /3 6.7, from which we estimate that in the 
chiral limit the chiral phase transition is at /3 = fi x = 6.7(1). The states with complex Wilson lines 
disorder into states with a real negative Wilson line close to this chiral transition. 



3. Nf = 3 simulations and results 

3.1 N, = 4 

We simulate lattice QCD with 3 flavours of colour-sextet quarks on 12 3 x 4 lattices, with 
quark masses m = 0.005, m = 0.01 and m = 0.02. j8s are chosen covering the range 5 < j8 < 7. 
Figure || shows the Wilson Lines and chiral condensates for a series of runs starting from an ordered 
start (gauge fields £/ = l)at/3=7. Note that the deconfinement and chiral-symmetry restoration 
transitions appear far apart. 

From the time evolution of the Wilson lines close to the deconfinement transition, our prelimi- 
nary estimates for the j3 values for this transition are fa{m = 0.02) = 5.295(5) and fa(m = 0.01) = 
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states with real positive Wilson Lines for Nf = 3 on 
a 12 3 x 6 lattice. 



5.285(5). (Pd(fn = 0.005) has yet to be determined.) We estimate that the chiral-symmetry restora- 
tion phase transition at m = is at j8 = [5 X = 6.0(1), from the peaks in the chiral susceptibilities. 
We are also performing a series of runs starting from a state where all gauge fields U = 1 except 
for the timelike fields on one timeslice which are set equal to diag(l, — 1, — 1) at j8 =7. This forces 
the system into a state with a negative Wilson line at large /3 s. From these runs we observe that 
states having Wilson lines with arguments ±2n/3 disorder into a state with a negative Wilson line 
above some /3 in the range 5.5 < j8 < 5.6. 



3.2 N t = 6 

We simulate lattice QCD with 3 flavours of colour-sextet quarks on 12 x 6 lattices, with 
quark masses m = 0.005, m = 0.01 and m = 0.02. j8s are chosen covering the range 5.3 < j8 < 7.0. 
Above the deconfinement transition the Wilson line exhibits a clear 3-state signal, so we bin the 
'data' appropriately. Figure || shows the Wilson Lines and chiral condensates for the states with real 
positive Wilson Lines for Nf = 3 on a 12 3 x 6 lattice, for all 3 masses. Again, the deconfinement 
and chiral-symmetry restoration transitions are well separated. 

By looking at histograms of the magnitude of Wilson Lines we estimate that the deconfinement 
transition for m = 0.02 is at fi d = 5.41(1). For m = 0.01 we estimate that f5 d = 5.395(5). We have 
yet to complete our simulations at m = 0.005 in the neighbourhood of the deconfinement transition. 
From the peaks of the chiral susceptibilities for each of the 3 masses, we estimate that the position 
of the m = chiral phase transition is at j8 = fi x = 6.3(1). There also exists a state with a negative 
Wilson Line at large j8. For some j3 in the range 6.1 < j8 < 6.2 this undergoes a transition to a 
state with its Wilson Line oriented in the direction of one of the complex cube roots of unity. At j3 
values just above the deconfinement transition we observe tunnelings between states with Wilson 
lines oriented in the directions of the 3 cube roots of unity. 
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4. Discussion and Conclusions 



We are simulating the thermodynamics of lattice QCD with 2 flavours of staggered colour- 
sextet quarks, whose chiral limit is a model of walking/conformal Technicolor. In addition we are 
simulating the 3-flavour theory, which is believed to have an infrared-stable fixed point, making it 
a conformal field theory at zero quark mass. 

We are extending our earlier Nf = 2 simulations at N t = 4, 6, to N t = 8. We again find widely 
separated deconfinement and chiral-symmetry restoration transitions. The /3 values for these tran- 
sitions are given in table [jj 
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Table 1: Nf = 2 deconfinement and chiral 
transitions. 



Table 2: Nf = 3 deconfinement and chiral 
transitions. 



Both j3rf and j8 x increase significantly with increasing N t . However, the fact that the increase 
in (5 X between N t = 6 and 8 is much smaller that that between N t =4- and 6 indicates that this is not 
simply due to the expected running of the bare (lattice) coupling with lattice spacing, controlled 
by asymptotic freedom. The most likely explanation is that for N t = 4 and 6, j8 = (5 X is in the 
strong-coupling domain. Here, once the condensate forms, the quarks play little part in the running 
of the coupling constant, which now runs as in the quenched theory. Between N t = 6 and 8 the 
theory is starting to emerge from the strong-coupling domain, and the effects of the fermions are 
starting to be felt. We will therefore need to simulate at even larger N t to see if the running of the 
coupling is now controlled by the asymptotically free g = fixed point or if a bulk chiral transition 
occurs before one reaches this fixed point. In the first case the continuum theory is confining and 
chirally broken, but with a slowly running - walking - coupling. In the second case the continuum 
theory is conformal. 

DeGrand, Shamir and Svetitsky have studied QCD with 2 sextet Wilson quarks They also 



find it difficult to distinguish between walking and conformal behaviour [JlOp. At finite temperature, 
they do not, however, see any separation between the deconfinement and chiral transitions JIT|]. The 
Lattice Higgs Collaboration 1 have been studying QCD with 2 sextet improved staggered quarks at 
zero temperatures. Their preliminary results are consistent with the walking scenario [|T2|]. 

We are also simulating the thermodynamics of lattice QCD with 3 flavours of staggered colour- 
sextet quarks. Since this theory is expected to be conformal, it is interesting to compare its be- 
haviour with that of the 2-flavour model. In our present simulations with N t = 4, 6, the results look 
very similar to the Nf = 2 theory except that the transitions are (as expected) moved to smaller 
j8s. We have tabulated these results in table ^. The deconfinement and chiral transitions are well 
separated and move to significantly higher j8 values as N t is increased from 4 to 6. Since this the- 
ory is expected to be conformal, the large increase in fi x is almost certainly quenched scaling in 
the strong-coupling regime. In fact, fi x at N t = 4 is close to for the quenched theory at N t = 8 



Z. Fodor, K. Holland, J. Kuti, D. Nogradi, C. Schroeder. 
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while that at N t = 6 is close to the quenched /3j at N t = 12 [ p"3| , [14], [L5p. Hence j3 x does appear to 
evolve with quenched dynamics. Here we would expect that, once we move to large enough N t , 
fi x will approach a fixed finite value representing a bulk transition separating the strongly coupled 
chirally-broken phase from the weak-coupling conformal phase. 

Above the deconfinement transition we see a 3-state signal, the remnant of the now-broken 
Z3 centre symmetry, for both Nf = 2 and Nf = 3. At an even higher j3, the 2 states with complex 
Wilson Lines disorder into a state with a negative Wilson Line. The existence of states with Wilson 
Lines having phases ±2n/2> and 71 in addition to those with phase is predicted by Machtey and 



Svetitsky and observed in their simulations with Wilson quarks [16]. 

As well as needing lattices with larger N t , we need to simulate at larger spatial volumes, 
especially for j8 in the neighbourhood of fi x and at low quark masses, to understand finite size 
effects. We are currently extending our Nf = 2, N t = 8 simulations with m = 0.005 to 24 3 x 8 
lattices. Nf = 2, N t = 12 simulations are planned, as are Nf = 3,N t = S simulations. These will 
concentrate on the vicinity of the chiral transition. 

Smaller quark masses are needed to be certain that we can access the chiral limit. Zero tem- 
perature simulations are needed to study hadron spectra, f K , interquark potential(s)..., to better 
determine whether these theories are walking or conformal. Such simulations will also allow us 
to study the flavour-symmetry ('taste') breaking from the use of staggered quarks, at relevant j8 
values. Direct measurement of the running of the coupling constant should be performed, along 
with measurement of the anomalous dimension of the chiral condensate. 
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